Visibility stacking in the quest for SNIa radio emission. 
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ABSTRACT 

We describe the process of stacking radio interferometry visibilities to form a deep com- 
posite image and its application to the observation of transient phenomena. We apply "vis- 
ibility stacking" to 46 archival Very Large Array observations of nearby type la supernovae 
(SNela). This new approach provides an upper limit on the SNIa ensemble peak radio lumi- 
nosity of 1.2 X 10^^ergs~^Hz^^ at 5 GHz, which is 5-10 times lower than previously measured. 
This luminosity implies an upper limit on the average companion stellar wind mass loss rate of 
1.3 X lO~^M0yr~-'^. This mass loss rate is consistent with the double degenerate scenario for 
SNela and rules out intermediate and high mass companions in the single degenerate scenario. 
In the era of time domain astronomy, techniques such as visibility stacking will be important in 
extracting the maximum amount of information from observations of populations of short lived 
events. 

Subject headings: stars: supernovae: general — techniques: interferometric 



1. Introduction 

Type la supernovae (SNela) have been impor- 
tant to astronomy since their identification as 
standard candles, and subsequent use in mea- 
sureing the accelerating expansion of the u niverse 



(|Riess et all Il998t iPerlmutter et all Il999l ). The 



current uncertainties in measuring cosmological 
parameters are domin ated by systemat ic rather 



than statistical errors (jNeill et al.l 120061 ). Hence, 



in order to better confine the cosmological pa- 
rameters it is important to understand the origin 
and evolution of SNela (jHillebrandt fc Niemever 
20001) . 



Recent work has been of an empirical nature, 
characterizing the differences betwe en SNIa light 
curve s as a function of environment (iHamuv et al 
2000 ). redshift, and metallicity (jTimmes et al 
20031 ). Another approach is to understand the 



physics of a SNIa explosion from a fundamen- 
tal level so that the shape and luminosity of the 
light curve, and its subsequent evolution, can be 
predicted and understood more completely. One 



of the key parameters affecting the light curve 
evolution is the amount of circumstellar material 
(GSM) that is present. 

As well as causing optical extinction, the ex- 
pected GSM provides a medium with which the 
supernova (SN) ejecta can interact and produce 
radio synchrotron emission. This process observed 
in type II and Ib/c SN. The spectral and tempo- 
ral distribution of this emission has been used to 
explore the envi ronment and histor y of the pro- 
genitor systems ( Montes et al.lll99"8l ). Despite the 
usefulness of this approach in characterizing the 
progenitor system, the radio detection of SNela 
has remained elusive. 

The mechanism by which SNIa white dwarf 
progenitors obtain a critical mass of ^ 1.4Mq is 
thought to be either accretion from, or merger 
with, a binary companion. In the single degener- 
ate (SD) scenario, a GO white dwarf accretes ma- 
terial from a companion star which is losing mass 
via Roche-lobe overf l ow or strong stellar wind 
(|Whelan k Ibenlll973tlNomotolll982h . In the dou- 
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ble degenerate (D P) scenario, two white dwarfs co- 
alesc e and ignite (|Webbinklll98lllben fc Tutukovl 
198J). In both cases the strength of the radio 
emission depends on the density of the CSM which 
is in turn dependent on the density of the progen- 
itor stellar wind. The shocked CSM can also pro- 
duce X-r ay emission which c an be used to probe 
the CSM. Ilmmler et alJ (|2006l ) claimed to have de- 
tected X-ray emission from the shocked CSM of 
th e SNIa 2005ke , howe ver this was not confirmed 
bylingheseEaD (|2007l) . 

Many projects have searched for radio emis- 
sion from SNela, but no detection has bee n made 
(|Eck et alJliooa Il995t IWeiler et al.lll989l l. The 
largest s et of radio observati ons of SNela is col- 
lected bv lPanagia et alJ (|2006l) . In the SD scenario 
the current lack of a radio detection has placed 
constraints on the companion mass loss rate, and 
restricted the mass of a post main-sequence com- 
panion to be relatively low. At even lower lu- 
minosity and mass loss rates, the solely-SD sce- 
nario becomes untenable, whilst a radio detection 
of a SNIa would dis-confirm the solely-DD sce- 
nario which predicts no such emission. Radio ob- 
servations of SNela can thus be used to determine 
whether or not all SNela result from only the SD 
or DD scenario, and can probe the ratio of SD to 
DD progenitors, in the case of a mixed population. 

More sensitive measurements can be obtained 
from deeper observations, or via a combination 
of currently available observations. Stacking is 
the technique of combining observations or images 
to form a composite, and is the two dimensional 
equivalent of forming a weighted average. Stack- 
ing can been used to extract the average proper- 
ties of a population of sources that are too faint 
to be de tected individually and has been used 
at 7- ray ( Aleksic et al .''201lV X-ray ( Laird et al 



2010h. optical (IWhite et al..,2007l) . infrared and ra- 
dio (|Garn fc Alexandeijl2009l )~ wavelengths. 

This letter is organized as follows: ^ outlines 
the process of stacking and introduces visibility 
stacking, ^J3] describes the data selection and re- 
duction process, §4]-[5] present and discuss the re- 
sults of the visibility stacking, and conclusions are 
drawn in ^ 



2. Visibility Stacking 

Traditionally stacking involves forming cali- 
brated images of a single source which are then 
weighted and combined to produce a more sen- 
sitive image. For N observations, each with an 
integration time of t, and Gaussian noise with an 
rms (T, a weighting scheme of w — will result 
in a combined image with rms noise of a/^/N, 
which is equivalent to a single observation of to- 
tal length T = N ■ t. Such a stacking scheme is 
commonly used to break long observations into 
multiple shorter ones, to mitigate against prob- 
lems such as radio interference, cosmic rays, and 
time-dependent calibration schemes, depending 
on the wavelength of observation. Non-Gaussian 
or correlated noise signals will result in a stacked 
image with a noise that is above the ideal a /^/N , 
and different weighting schemes are required to 
account for such problems. 

Imaging radio synthesis observations involves 
a conversion between the visibility and image 
domains. The visibility sampling function (or 
{u, v) coverage) of a radio interferometer is in- 
complete and generally non-uniform. Incomplete 
(w, v) coverage will cause the reconstructed im- 
age to deviate from the true radio sky, with 
some spatial scales being under represented or 
absent. This problem is well known and moti- 
yates the design a nd layout of radio observatories 
( Thomps"onl Il999) including the Very Large Ar- 
rav fVLA. [Thompson et al.lll98"ol). Australia Tele - 
scope Compact Array fATCA. lFrater et al.lll992l) 
the Low Frequency Array (LOFAR. iRottgerind 
20031 ) and the Australia n Square Kilometer Ar- 
ray Pathfinder (ASKAP. Ijohnston et aLll2009t ). A 
common method for combating incomplete (u, v) 
coverage is to combine observations that have 
complimentary coverages. For fixed arrays 

such as LOFAR and ASKAP this can be achieved 
with earth rotation synthesis, while movable ar- 
rays such as the VLA and the ATCA can also 
use multiple configurations. This commonly used 
process is the trivial case of stacking observations 
of a single object that are temporally separated. 

Stacking of radio images of a population of 
sources can be done in the same manner as 
stacki ng at other wavelengths (in the ima^e do- 
main, 



Garn fc Alexandeijl2009l : iHodge et al.lboOSt 

White et al.ll2007h . however incomplete (u, v) cov- 
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erage can cause problems. Observations that are 
obtained from different configurations of a tele- 
scope, will have a different (u, v) coverage and 
thus be sensitive to a different range of spatial 
scales and have different resolutions. An isolated 
point source can be easily detected at all spatial 
scales so stacking in the image domain will still 
recover the source. However, the noise statistics 
of the image is not easily understood. In more 
complicated fields, the distribution of sources and 
their morphologies become confused with spatial 
sampling effects and confident detections are dif- 
ficult. 

It is possible to benefit from improved (u, v) 
coverage as well as the lower noise gained from 
a longer integration time by using a technique 
we term visibility stacking. Visibility stacking is 
achieved by combining visibilities from observa- 
tions of a population of sources before the image 
is created. 

When stacking in the image domain it is im- 
portant to ensure that the input images are free 
of confusing sources, and the target source posi- 
tions are aligned. The same is true for visibility 
stacking, but these corrections need to be applied 
to the visibility data. Removing confusing sources 
can be achieved by subtracting appropriate mod- 
els from the visibility data. Aligning the target 
source positions can be achieved by adjusting the 
phases of the visibilities. 

3. Target selection 

Observations of 27 nearby (z < 0.1) S Nela have 
been carried out bv lPanagia et al. I (|2006l) using the 
VLA at wavelengths between 0.7 and 20 cm from 
1981 to 2003. This data set has a large number 
of SNIa observations from a single instrument and 
a consistent set of calibrators making it ideal for 
forming a stacked image. 

The data set is not completely homogeneous, 
with the VLA configuration being restricted by the 
observing season. The archival observations cover 
all standard configurations of the VLA (A, B, C, 
D, AnB, BnC, CnD), and are spread across mul- 
tiple frequencies. In order to obtain the most sen- 
sitive stacked image, the 6 cm observations were 
used as they constituted 75 of the 148 observa- 
tions. 



3.1. Data reduction 

The raw data were obtained from the VLA 
arch ivj^, converted into uvfits format wi th AIPS 
JCreisen 2003), and reduced in MIRIAD (jSault et al, 
Il995[ ). Flagging was based on manual inspec- 



tion. I n keeping with the analysis of lPanagia et al 



(120061 ). 3C286 was used as a primary calibrator for 
each of the observations. The flux of the primary 
calibrator varied by less than 1% over the course 
of the observing period. 

A number o f obse rvations that were listed by 
Panagia et al. ( 2006f ) were not used in the final 
stacked image due to problems with: data access 
(4 could not be located in the VLA archive) , data 
quality (6 had poor or missing calibrator observa- 
tions, 4 had incorrect gains), or complicated field 
sources (15 observations, 8 of which include Cen- 
taurus A). The total number of usable observa- 
tions was thus reduced to 46, and are listed in 
Table [H 

In many cases the field of view of the VLA 
included background radio sources that created 
side-lobes and artefacts at the expected SNIa po- 
sition which needed to be removed. A model of 
the sources was removed from affected visibilities. 
Column 4 of Table [1] lists the rms noise of indi- 
vidual observations when naturally weighted and 
with background sources removed. 

Visibility stacking was achieved in MIRIAD by 
adjusting the phase center of each image to be at 
the coordinates of the target SN, and imaging with 
the task invert. Many of the observations con- 
tained two 50 MHz bands, 4.885 and 4.835 GHz 
(and one observation at 4.889 GHz), each of which 
were calibrated separately. A total of 79 such 
bands were used in the creation of the stacked im- 
ages. 

3.2. Interpretation 

Images created via visibility stacking require 
some interpretation. In order to retain the great- 
est image sensitivity, the visibilities were naturally 
weighted when imaging. The flux in the image can 
thus be interpreted as being from a source at an 
effective distance do, which is a weighted average 
of the input source distances. The weighting func- 
tion for do is then the same as for the visibility 



^ https : / / archive, nrao.cdu /archive /advquery. j sp 
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Tabic 1: Archival VLA observations used in the stacking analysis. 



Name 


Obs Date 


Age 


fecmA'Jy 


dist 


Name 


Obs Date 


Age 




dist 


SN 


Y-M-D 


days 


beam~^ 


Mpc 


SN 


Y-M-D 


days 


beam~^ 


Mpc 


1981B 


1981-Mar-ll 


18 


82 


16.6 


19860 


1987-Fcb-12 


72 


57 


28 


II 


1981-Apr-9 


46 


105 


II 


II 


1987-Apr-ll 


130 


68 




II 


1982-Jun-25 


489 


31 


II 


II 


1987-May-24 


173 


67 




II 


1982-Oct-2 


587 


50 


II 


II 


1987-Aug-28 


269 


52 




1982E 


1985-Dec-29 


1417 


44 


23.1 


II 


1989-Jul-17 


958 


90 




1983G 


1983-May-27 


72 


131 


17.8 


1987D 


1987-May-15 


46 


108 


30 


1984A 


1984-Mar-5 


430 


109 


17.4 


II 


1987-Jun-4 


66 


124 




1985A 


1985-Feb-l 


49 


93 


26.8 


II 


1987-Jun-21 


83 


32 




II 


1985-Feb-17 


65 


65 


II 


II 


1987-Sep-18 


172 


62 




II 


1986-Jun-15 


550 


71 


II 


II 


1988-May-29 


426 


174 




II 


1987-Oct-23 


1045 


104 


II 


1987N 


1987-Dec-20 


37 


91 


37.0 


1985B 


1985-Feb-22 


70 


173 


28.0 


II 


1988-Feb-l 


76 


84 




II 


1985-Mar-18 


94 


50 


II 


II 


1988-Mar-31 


135 


55 




II 


1985-Sep-15 


275 


74 


II 


II 


1988-Aug-22 


279 


64 




II 


1986-Apr-30 


502 


81 


II 


1989B 


1989-Feb-2 


10 


38 


11.1 


II 


1987-Scp-18 


1008 


56 


II 


II 


1989-Feb-3 


11 


27 




1986A 


1986-Feb-7 


18 


126 


46.1 


II 


1989-Mar-6 


42 


39 




II 


1986-Feb-25 


36 


89 


II 


II 


1989-Mar-27 


63 


57 




II 


1986-Mar-16 


55 


38 


II 


II 


1989-Apr-6 


73 


64 




II 


1986-Apr-3 


73 


51 


II 


1989M 


1990-Feb-13 


271 


64 


17.4 


II 


1986-Jun-15 


146 


67 


II 


II 


1990-May-29 


376 


65 




II 


1987-Apr-l 


436 


68 


II 


1992A 


1993-Feb-5 


400 


53 


24.0 












1994D 


1994-May-4 


57 


35 


14 












199Sbu 


1999-Jaii-7 


24.J 


(jT 


11.8 
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Table 2: Parameters for the stacked images. Radio luminosity and inferred mass loss rates are 3cr upper 
limits. 







Image Stack 






Early 


Late 


Ah 


do (Mpc) 


14.3 


19.0 


15.7 


Effective Age (Days) 


51.9 


478 


251 


Effective Integration 
at 50 MHz b/w (hr) 


22.5 


19.6 


42.1 


la (/xJy beam^^) 


16.7 


19.2 


13.3 


Synthesized beam (arcsec^) 


21 X 18 


22 X 15 


21 X 17 


Lo(ergs-iHz-i) 


1.2 X 10^^ 


2.5 X 10^^ 


1.2 X 10^5 


M(Moyr-i) 


1.3 X 10-^ 


1.5 X 10-6 


5.4 X 10-^ 



data: 



Ni 



N2 



N„ 
"5E 



Ni+N2 + --- + N„ 



(1) 



where Ni is the number of visibilities in i*'^ ob- 
servation, and di is the distance to the correspond- 
ing source. The function is of the form Ni/df since 
doubling the integration time doubles Ni and in- 
creases the sensitivity by which allows a source 
of the same luminosity to be seen at a distance 
times further away. 

Using Equation ([T]), it is possible to calculate 
the effective distance probed by a stacked image, 
and convert the flux density of a detection or up- 
per limit into an ensemble luminosity or luminos- 
ity limit. Distances listed in Table [T] are from 



Panagia et al.l (|2006l ). and were used to calculate 
the effective distance do for the stacked image. 

A similar weighting function was applied to the 
age of a SN at the time of its observation to obtain 
an effective age that corresponds to the luminosity 
calculated in the above-mentioned manner. The 
relevant weighting function is: 



Age = 



Agci • A^i + Agcj • N2 



Age„ • N„ 



Ni+N2 + --- + N„ 



(2) 



Results 



The 46 usable observations were combined to 
make three stacked images of the SNe. Two im- 
ages were formed by selecting only observations 
with an age of either less than 100 days (Early) 
or greater than 100 days (Late), respectively. A 



third image was created using all of the available 
observations (All). In all three images no emis- 
sion was seen above 3(j at the location of the SNe. 
The representative distance, age and correspond- 
ing luminosity of each of these images are listed 
in Table [21 and are plotted in Figure [1] along with 
the upper limits obtained from the individual ob- 
servations. 

The stacked image of all observations represents 
a 50 MHz band width equivalent integration time 
of 42.1 hours, longer than any previous radio ob- 
servation of any individual SNIa. 

The final stacked images were not completely 
without structure and some faint sources were 
detected away from the expected SN location. 
Stacked images were created for each of the host 
galaxies, and faint sources that were not able to 
be seen or removed from the short observations 
were detectable in the more sensitive stacked im- 
age. These sources were found to be HII regions 
within, and the core of, NGC 362 7 and were pre- 
viously note d bylEck et al.l (|2002l). 

Following lPanagia et al. (2006, 2011 ). we model 
the luminosity of the SNla as: 



T = A I 

1026ergs-iHz-^ I 



/ M/lO^^Moyr- 



1.65 



5GHz 



Wwind/lOkms 

-1.5 



1 day 



(3) 

where t is the time since explosion, or age of the 
SN. The parameter A is a dimensionless constant. 
On the assumption that the light curves of SNela 
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Fig. 1. — The Sa noise of the individual and stacked images as a function of SN age and distance. Three 
stacked images were created - Early, All and Late. In each stacked image the noise level is fainter than in 
the individual observations. Left: The noise in each image as a function of the age of the SNe at the time of 
observation. Right: The noise in each image as a function of SN distance. Lines have been drawn to show 
the 3cr upper limit for different 5 GHz luminosities in ergs~^Hz^^. 



are similar to that of SNelb/c, we take A = 1285 ± 



245 as calculated by IPanagia et al] (|2Q06l ). The 
value of A can vary by a factor of 10 depending 
on which set of SNe are used in the derivation, 
however the dependence of Lj^ on M/wwind results 
in an uncertainty in M of less than a factor of 2. 
A wind velocity Wwind = 10km s^^ was assumed, 
consistent with a post main-sequence companion 
in the SD scenario. This model assumes that the 
CSM has a radial density profile p oc r^^. 

The optical depth tcsm of the CSM external 
to the blast wave depends on absorption mech- 
anism (free-free or synchrotron self absorption), 
as well as the structure of the CSM (clumpy or 
uniform). Constant among all considerations of 
optical depth is that the luminosity peaks when 
the optical depth is unity. With the assumption 
'''CSM = 1, equation ([3|) was used to estimate an 
upper limit for the mass loss rate M, as listed in 
Table [I 

5. Discussion 

Combining data from different observations is a 
common technique for making measurements for a 
population of sources, however when dealing with 
upper limits care needs to be taken. A combined 



upper limit on the radio luminosity of SNela stud- 
ied here can be obtained by creating a weighted av- 
erage of the individual upper limits without stack- 
ing. Whilst this will result in a combined upper 
limit that is consistent with that obtained from 
a stacked image, it does not allow for the detec- 
tion of sources. Stacking images permits the de- 
tection of source which would otherwise be below 
the detection limit, whilst combining upper limits 
implicitly assumes that no such detection will be 
made. 

The stacked 5 GHz images provide upper lim- 
its to the ensemble radio flux of SNela which 
are fainter than the individual observations. The 
stacked images of Early and All observations both 
have a Scr upper limit on the 5 GHz luminosity of 
1.2 X lO^^ergs^^Hz^^ which is the same as that 
of the best observation of SN1989B (1989-Feb-03). 
Whilst the stacked images are more sensitive than 
the individual observations (see Fig [T] left), the 
fact that (X (jcfi means that the flux/distance 
combination of the single observation and stacked 
image both correspond to the same luminosity 
limit (1.2 X IG^Sergs-^Hz"^ Fig [1] right). Rear- 
ranging Equation |3] with L^, cx a (P , we find that 
M cx a'^-^^d^-^t'^'^^ , so that the limit on mass loss 
rate is more strongly dependent on the distance of 
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the SN and age of observation than its sensitivity. 
Since the observation of SN1989B was made when 
it was only 10 days old, and the source is at only 
11.1 Mpc, it provides a lower limit on the com- 
panion mass loss rate (3.1 x lO^^MQyr"^) than 
the older and more distant Early stacked image 
(1.3x 10~^MQyr~-^) even though the Early stacked 
image is more sensitive. 

However, while the SN1989B observation pro- 
vides an upper limit on SN luminosity which is 
the same as Early and All stacked images, it is 
only applicable to SN1989B and is not represen- 
tative of the underlying population. If SNela are 
a completely homogeneous population of sources 
then this limit may be applied to all SNela, how- 
ever the opt ic al lum inosit y differences s e en b; 



Hamuv et ahl (|2000l ) and Iximmes et all (fm 



suggest that this is not the case. By creating 
a stacked image of the SNela we are able to 
claim a luminosity limit which is applicable to 
all SNe within the sample, and be more confi- 
dent that it is applicable to SNela as a whole. 
The 3tT lu minosity limits and mass loss rates re- 
ported by iPanagia et al. (i2006jl have a median 
of 1.2 X 1026ergs-iHz-^ and 6.2 x lO-^Meyr"! 
respectively. The ensemble luminosity limits for 
the Early, Late and All stacked images are 5 — 10 
times lower these earlier limits, and the ensemble 
mass loss rate constraints of the Early and All 
stacked images are up to 5 times lower. 

If the solely-SD scenario is correct then the 
mass transfer from the companion to the WD 
must result in some residual mass that is not ac- 
creted, our observations limit this contribution to 
the CSM to a rate of less than 1.3 x 10~'^MQyr"^ 
Wind accretion onto a binary companion is typi- 
cally ^ 10% efficient (lYungelson et al.lll995l ). im- 
plying a companion mass loss rate of < 1.4 x 
lO~^M0yr~^ ruling out red giants and medium to 
high mass (2 — IOMq) post main-sequence stars as 
the possible binary companion. A main-sequence 
companion, having winds 3> 10km s'^, could re- 
main undetected even with much higher mass loss 
rates. Accretion via Roche-lobe overflow (RLOF) 
is mor e efficient than wind a ccretion and the mod- 
els of Hachisu et al. (Il999h predict that SNela 
formed from RLOF of a 1.3— 2M0 red giant on to a 



^ Their 2a upper limits have been converted to 3cr limits for 
comparison. 



CO white dwarf are consistent with the measured 
mass loss rate. The solely-SD scenario (with a red 
giant or post main-sequence companion > 2Mq) 
is thus not able to account for the properties of the 
observed SNIa population. A detailed analysis of 
the formation of binary systems with large mass 
ratios (5 — SM© primary, < 2Mq secondary) is 
beyond the scope of this paper, but such systems 
may be relatively rare because of the amount of 
mass transfer from the primary to the companion 
during the final stages of the its evolution. 

6. Conclusions 

By employing the new technique of visibility 
stacking we were able to create more sensitive im- 
ages of SNela. The corresponding 3ct limits on the 
ensemble radio luminosity are considerably lowe r 
than most of the limits of IPanagia et al.l ( 2006 ). 
Since the calculated mass loss rate depends on the 
distance and age of the SN as well as the sensi- 
tivity of the observation, we were unable to ob- 
tain an ensemble limit on M that was lower than 
the individual observation of SN1989B on Feb 3'"'' 
1989. However, we were able to place an overall 
stringent limit on the ensemble 5 GHz radio lumi- 
nosity and mass loss rate for the SNIa population 
of 1.2 X lO^^ergs-^Hz"^ and 1.3 x lO-'^Mgyr"! 
respectively. By creating visibility stacked images 
of SNela we conclude that the solely single de- 
generate scenario (with a red giant or post main- 
sequence companion > 2Mq) cannot account for 
the observed SNIa population. 

Visibility stacking is a powerful tool for analyz- 
ing radio observations of transient sources. The 
use of visibility stacking can be extended to the 
radio properties of any sample of radio quiet or 
radio faint sources, where image domain stacking 
has previously been used. The increased sensitiv- 
ity can increase the confidence with which detec- 
tions and upper limits are made. This technique 
needs to be applied carefully in the analysis of 
transient sources when such phenomena may re- 
sult from more than one progenitor scenario. 

In the era of time domain astronomy there will 
be a focus on obtaining more sensitive measure- 
ments on populations of sources that are transient. 
Such sources do not afford the luxury of repeated 
observations and thus techniques such as visibility 
stacking will be important in investigating source 
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populations. 

In studies which hope to detect the radio emis- 
sion of SNela by increasing; the sensiti vity of the 
observations (eg, IChomiuk et al. 2011 ) , the abil- 
ity to form visibility stacked images will further 
increase the sensitivity and applicability of the ob- 
servations. 
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